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We first develop a method to calculate a complete set of mode functions which describe the 
• quantum fluctuations generated in one-bubble open inflation models. We consider two classes of 

f"^ | models. One is a single scalar field model proposed by Bucher, Goldhaber and Turok and by us as an 

example of the open inflation scinario, and the other is a two-field model such as the "supernatural" 
inflation proposed by Linde and Mezhlumian. In both cases we assume the difference in the vacuum 
energy density between inside and outside the bubble is negligible. There are two kinds of mode 
functions. One kind has usual continuous spectrum and the other has discrete spectrum with 
characteristic wavelengths exceeding the spatial curvature scale. The latter can be further devided 
into two classes in terms of its origin. One is called the de Sitter super-curvature mode, which arises 
O |' due to the global spacetime structure of de Sitter space, and the other is due to fluctuations of the 

bubble wall. We calculate the spectrum of quantum fluctuations in these models and evaluate the 
resulting large angular scale CMB anisotropies. We find there are ranges of model parameters that 
are consistent with observed CMB anisotropies. 
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I. INTRODUCTION 

Motivated by observations that we may live in a low density universe [P}t3|, several authors have considered a 
possible scenario which realizes an open universe (f^o < 1) in the context of inflationary cosmology [0-0- In contrast 
to the standard inflation models which predict a spatially flat universe (Oo = 1), this scenario predicts an universe 
with spatially negative curvature. In this scenario, the nucleation of a vacuum bubble plays an essential role. In 
general, the bubble nucleation process is described by the bounce solution with 0(4)-symmetry, which is a non-trivial 
classical solution of the field equation in Euclidean spacetime ||J|. Then the expanding bubble after nucleation is 
described by the classical solution obtained by analytic continuation of the bounce solution to Lorentzian spacetime. 
Owing to the 0(4)-symmetry of the bounce solution, the expanding bubble has the 0(3,l)-symmetry. This implies 
that the system is homogeneous and isotropic on the hyperbolic time slicing inside the bubble and that the creation 
of one bubble can be regarded as the creation of an open universe. 

To obtain a realistic model of an open universe, the flatness f^o ~ 1 and the homogeneity and isotropy of the universe 
should be realized inside the bubble. This requirement can be satisfied by assuming, for example, a scalar field with 
a potential such as in the new inflation scenario but with a high potential barrier before the slow rolling inflationary 
phase In addition, the high potential barrier keeps the bubble collision rate small, hence the homogeneity and 

isotropy of the one-bubble universe is not disturbed by the other nucleated bubbles. Then the second inflation inside 
the bubble inflates the universe and explains the flatness £Iq ~ 1. Linde and Mezhlumian has proposed another 
(perhaps more natural) model of open inflation by introducing two scalar fields ^jj7|, but the essential feature of 

scenario is the same as the former. 

The interesting problem is the origin of fluctuations in this open inflationary scenario. Following the usual picture 
that quantum fluctuations of the scalar field generate the initial density perturbations, we need to investigate the 



quantum state of the scalar field after bubble nucleation. According to previous investigations, the open inflation 
scenario shows interesting varieties of the fluctuations. The bubble nucleation process can excite fluctuations of the 
scalar field and may increase the power of density fluctuations on scale of spatial curvature [|lO|]f l| ]. Also it has been 
shown that peculiar discrete modes of fluctuations on super-curvature scale may exist [fl2f and contribute to cosmic 
microwave background (CMB) anisotropics in an open universe f^,[l3|j. Very recently, generation of another type of 
super-curvature perturbations which originate from the bubble wall perturbations has been discussed pl|,p^ , ^5[ . 

In this paper, we develop a method to calculate these varieties of fluctuations in the open inflationary universe 
and evaluate the power spectrum of the resulting CMB anisotropics. In order to perform a detailed analysis of the 
power specturm, it is necessary to specify the model to some extent. Here we consider two classes of models. One 
is a single scalar field model (Model A) proposed by Bucher, Goldhaber and Turok [Ej and by ourselves ||. The 
other is a two-field model (Model B) in which the false vacuum decay and inflation inside the bubble are governed by 
two different scalar fields, such as the supernatural inflation proposed by Linde and Mezhlumian For simplicity, 
we assume the difference in the vacuum energy density between inside and outside the bubble is negligible in both 
models. Our method is based on the formalism we developed previously for computing the mode functions during 
and after bubble nucleation fl6|, 17|Jl^| . 

This paper is organized as 10 hows. In section 2 we describe our method to calculate quantum fluctuations of the 
inflaton field. We assume that the reaction to the geometry can be neglected, i.e., the nucleation occurs in the de 
Sitter spacetime background. Then, in section 3, we investigate the evolution of fluctuations inside the bubble in an 
open inflationary stage and derive the formulas to calculate the initial perturbation spectrum of the open universe and 
the resulting CMB anisotropics. In section 4 we evaluate the CMB anisotropics in two simple models and examine 
their viability. Section 5 is devoted to conclusion. We adopts the units c = K = 1, and the bar denotes the complex 
conjugate. 



II. FORMALISM 



In this section, we describe our formalism to investigate quantum field fluctuations inside a bubble. Pioneering 
work on the quantum state of a nucleating bubble was done by Rubakov |l8| , Kandrup (T^| , Vachaspati and Vilenkin 
p0|. Recently, we have developed a formalism to investigate the quantum state of a scalar field after false vacuum 
decay based on the WKB wave function of a multi-dimensional tunneling system p6| , pl| . The formalism has been 
applied to the bubble nucleation which occurs on the Minkowski background, and the spectrum of field fluctuations 
after the decay has been studied ]lO|,[Ll]]. The basic formalism has been extended to the case in which the bubble 
nucleation occurs on the de Sitter spacetime |17|. There the effect of the non-trivial geometry of the instanton with 
gravity, i.e., the Coleman-De Luccia instanton||] was taken into account. At that time, however, the appropriate set 
of the mode functions to describe the initial vacuum state, which is expected to be the Bunch-Davies vacuum due to 
the sufficient inflation before the tunneling, had not been known. Thus there was a technical difficulty in applying 
our formalism. Recently, we have succeeded in describing the Bunch-Davies vacuum state in the spatially open chart 
| p2[ . Thus this problem has been overcome. Combining the results in these two papers |l^,0, we now have a tool to 
bundle the quantum fluctuations after tunneling which includes the effect of the geometry ofthe Coleman-De Luccia 
instanton similar to the case in the Minkowski background fllCfl . 

In the present paper, we consider simple open inflation models in which the geometry of the Coleman-De Luccia 
instanton can be approximated by the pure de Sitter spacetime. That is we consider the case in which the potential 
energy difference between the two vacua of the tunneling filed is small. We consider the action 

S = I (A9^d^d u n - V(a)) ^-gd A x + / (-\g^d^d v 4> - U(a, ft) V^g^x, (2.1) 



2" -<*-~"- v 7 V a ~ ~ ' J \ 2 

where the potential V(a) is assumed to have the form as shown in Figs. I (a), (b) to realize the false vacuum decay. 
As we mentioned above, we take account of the effect of gravity only as a curved background which is assumed to 
be de Sitter space. We denote the value of a at the false vacuum by of- The field 4> is the inflaton in the nucleated 
bubble. We devide them into two parts as <p = 4>b + f- 4>b is the seimi-classical, spatially homogeneous part of the 
field and <p> represents the quantum fluctuations around this background. We neglect the quantum fluctuations of 
a and denote its 0(4)-symmetric background by <jb- Due to the 0(4)-symmetry of the instanton, ub is a function 
of only one coordinate, say r. In addition, we assume 4>b also respects the 0(4)-symmetry. Thus a model we have 
in mind is a two-field one, with <Pb being constant at <jb = °\f- Note that we can consider a single field model by 
identifying (p with the quantum fluctuations of a and letting U(a, ft — V" (<j)4> 2 /2. 
With these assumptions, the action for tp reduces to 



S, J [ -lg»»d»<pd v <p- ±M 2 . r)r ) \ —./',-. (2.2) 



where 
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^ 2 (*t) = S(^(t),<Mt)). (2.3) 



A2 

<2/ 



As noted above, in the case of a single field model, 4>b and (r) in the following discussion should be replaced by 
d 2 V 

gb and — 2 (<tb(t)) , respectively. 

Our formalism is summarized in the following. It is based on the WKB wave function which describes the tunneling 
decay in a multi-dimensional system . The construction of the wave function was originally developed by Gervais 

and Sakita p2]] . The wave functional is written in terms of the bounce solution as, which gives the semi-classical 
picture of the nucleation and expanding bubble, and a set of mode functions which describe the quantum fluctuations 
(p on the bounce solution |l6| . The boundary condition for these mode functions is determined from the fact that 
the wave function describes the false vacuum state at <jb — &f- It then turned out that the procedure to obtain 
the appropriate set of mode functions is equivalent to finding out a complete set of them which are regular in one 
hemisphere of the Euclideanized de Sitter space fl7j , |l2|| . That is, to find a complete set of mode functions v k which 
obey 



«k(i,x)=0, (2.4) 



in the Lorentzian region, and which are regular on the Imt < hemisphere. Here, it is important to note that "a 
complete set" means a set of all modes which are properly normalized by the Klein-Gordon norms on a Cauchy surface 
S of spacetime, 

(i>k,fv) := ~i I dT,^" {v k d„v k > - (<9„u k )vk'} 

= <W. (2.5) 

Once they are obtained, the quantum fluctuations of the field are described by the "vacuum state", such that 
flkl'I') = for any k where the fluctuating field is expressed as 



k 



v k a k + v k al), (2.6) 



in the Heisenberg representation, with a k and being the annihilation and creation operators, respectively. Thus 
the mode functions v k play the role of the positive frequency functions. 

To write down the equation for the mode functions, we introduce the coordinates in the de Sitter spacetime following 
Rcf. jl2). The four-dimensional Euclidean de Sitter space is a four-sphere. The metric is represented as 

ds% = H- 2 dr 2 + a 2 E (T){dp 2 + sin 2 p dti 2 )], (2.7) 

where — 7r/2 < r < n/2, < p < tt and cle(t) = H^ 1 cost. The mode functions are required to be regular on the 
< p < 7r/2 hemisphere. 

The coordinate systems in the Lorentzian de Sitter space are obtained by analytic continuation as 

t R = i(T-ir/2) (t R >0), 
rn = ip (tr > 0), 
t c = r (vr/2 > t c > -7r/2), 
Tc = i(p — 7r /2) (oo > rc > — oo), 
t L = *(-t - tt/2) (t L > 0), 

r L =*P (r L >0). (2.8) 

we find that each set of these coordinates covers three distinct parts of the Lorentzian de Sitter spacetime, which we 
call the regions R, C, and L (see Fig. 2). The metrics in these three regions are given by 

ds 2 R = -H~ 2 dt 2 R + a 2 (t R ){dr 2 R + sinh 2 r R dQ, 2 ), 
ds 2 c = H~ 2 dt 2 c + a 2 E (t c ){-dr 2 c + cosh 2 r c dQ 2 ), 

ds 2 L = -H- 2 dt 2 L + a 2 (t L ){dr 2 L + sinh 2 r L dfl 2 ), (2.9) 

respectively, where a(t) = iJ _1 sinht. We assign the region L to be in the false vacuum sea and the region R to 
describe the open universe inside the bubble. 
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The equation for the bounce in the Euclidean region is given by 



H^d_ 

a E dr 



3 da B \ 



V'(a B )=0. 



(2.10) 



The equation in the Lorentzian region is given by the analytic continuation of the coordinates as specified in Eqs. fl2.S| ). 
The fluctuation field if obeys, in the Euclidean region, 



1 



d 



d 



1 



where 



sin 2 p dp 



sin p 



dp 



M 2 {t) 



T 2 

• 2 1 

sm p 



v = o, 



(2.11) 



(2.12) 



and Ln is the Laplacian on a unit two sphere. 



Setting (p — ob(t) 1 Xp( T )fpi(p)Yi m (Sl), the equations for Xp and f p i become, respectively, 



2 1 



(l-x 2 



dx 2 



- 2x 



dx 



l-x 2 



M 
H 2 



Xp = 0, 



where x = sinr, and 



d 



sin p dp 



sin p 



d_\ , 1(1 + 1) 
dp 



sin 2 p 



(P 2 



0. 



(2.13) 



(2.14) 



The requirement that the mode functions be regular at p = rj, = t\r = fixes the form of f p i to be E^ ] 



r(ip + ; + i) 



p. 



-(-1/2 



r(tp + i) vimh7 ip " 1/2 



(cosh r) 



2 r(-ip+l) 
7rr(-ip + ? + 1) 



sinh' r 



smpr 



d(coshr)' \sinh 



(2.15) 



where r = ip in the Euclidean region. The gamma function factor in front is attached so that the functions 
Ypi m (r, CI) := f p i(r)Yi m (CY) become properly normalized harmonics on the hyperbolic slices in the region R or L: 



f°° 2 f 

/ dr sinh r / 



drsinh r / dflY p i m (r,Cl)Ypn> m t(r, SI) = 6(p —p')8u>6 v 



(2.16) 



As we have noted before, a complete set of mode functions should be properly normalized on a Cauchy surface. 
Therefore we consider the mode functions in the region C, since it is the region in which complete Cauchy surfaces 
exist. It should be noted that f p i play the role of the positive frequency functions there. Then the remaining tasks 
are (1) to find the possible spectrum of p that may give finite Klein-Gordon norms and (2) to properly normalize the 
mode functions if they are normalizable. 

To accomplish the first task, we now consider the equation for % p , (2.13). In the region C, the equation for \p 
remains the same because tc — t. Introducing the conformal coordinate £ by the relation x = sinr = tanh£, the 
metric in the region C becomes 



ds 2 = a|(£) (-dr 2 c + d£ 2 + cosh 2 r c dSl 2 ) , 
where = (i? cosh£) _1 and Eq. ( 2.13| ) is rewritten as 



di 



X P =P Xp, 



where 



u(0 



M 2 /H 2 - 2 
cosh 2 £ 



(2.17) 



(2.18) 



(2.19) 
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This resembles the Schrodinger equation with the potential £/(£) and the energy E = p 2 . Thus finding possible values 
of p is equivalent to solving the eigenvalue problem of the above equation. Since L/(£) — > at £ — > ±00, the spectrum 
of modes is continuous for p 2 > 0. On the other hand, if A4 2 < 2H 2 for some interval of £, the potential has a valley 
and some discrete modes which correspond to bound states may appear for p 2 = —A 2 < 0. In fact, it has been shown 
that there exists a discrete eigenvalue of A 2 in two different simple cases. One is a mode which appears when we 
consider a scalar field with a constant small mass (< \/2H 2 ) on the de Sitter spacetime p2[ . Its origin is intrinsic 
to the spacetime structure of the de Sitter universe. As long as the volume with A4 2 < 2H 2 is large enough, an 
analogous mode is expected to exist in general. Let us call it the de Sitter super-curvature mode.^J The other is the 
mode that appears when we cosider the fluctuations of the tunneling field itself pj. It describes the fluctuation of 
the wall. Since it turns out that the analysis of the wall fluctuation mode can be simplified compared with the other 
discrete modes, below we consider the continuous modes, the de Sitter super-curvature mode, and the wall fluctuation 
mode separately in sequence. For simlicity, we assume that M. is constant in the regions R and L and set 



M 2 = {™1 ln f> (2.20) 
1 AF m L . 



in the following discussion. 



A. Continuous modes 



We first consider the continuous spectrum mode functions of p 2 > 0. For these modes, it has been shown that the 
Klein-Gordon norms can be evaluated by the sum of those on the hyperbolic slices on the regions R and L [jl2j . Since 
in the regions R and L, M 2 is supposed to be constant, \p is given by an associated Legendre function P* P _ 1 / 2 ( Z ) or 
P~l^/ 2 ( z ) with v = y9/4 — A4 2 / H 2 . To construct the normalized mode functions, we begin with the following mode 
functions in the regions R and L: 



where 



'£ = ^xHtfl)Wr H ,n), 

xl L) (t L )Y plm (r L ,n), 



plm 



a(t L 



(2.21) 



■ cosh tR , 
cosh it , 



(2.22) 



with v' = v/9/4 — m 2 1 H 2 — 1/2 and p! — ^/9/A~ M 2 /H 2 — 1/2. The norms of these functions are give by 



oXp_-(R) 
dz R Xp 

,1 ( L ) 
dZL p 




— sinh 7rp, 

7T 



• sinh 7rp. 



(2.23) 



If we analytically continue Xp to the region L by solving Eq. (2.18) in the region C, xi> wm be expressed in terms 



(R) 



of a linear combination of Xv an d X^—p- Hence we consider the mode function, 

Xp = Xp ^ = a pX—p "I" PpXp ^ ■ 



(2.24) 



From the property of the analytic continuation E2] and the reality of the Eq. (2.18), there exists a symmetry, 



*The spatial harmonics behave as Y p i m (r) oc e~ r for p 2 > . As r — 1 corresponds to the scale of spatial curvature, fluctuations 
described by the harmonics with p 2 > represents those which decay exponentially on scales larger than the curvature scale. 
In contrast, harmonics with p 2 = —A 2 < behave as Yju m oc e^ A_1 ' r , hence describe fluctuations over scales larger than the 
curvature scale fcd,p4j. The name "super-curvature" originates from this fact. 
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Ap A — p i 

Ap ° A— p 5 



in the region C. Then we obtain 
Further, defining the Wronskian by 

ctp and P p are expressed as 



cu = e- 27Tp a 



o ) Pp — P—p- 



.. dv du 



U L) J P R) )) a _(U R \x { - L p ) )) 



((Xp ,X-p)) 



, pp 



U L \x w P )) 



and we have 



((xp,x P )) = ((xl R) ,Xp R) )) = (^ P K\ 2 ~ Wv\ 2 ) {{x { p L \x P L) )) 



Also from the analytic continuation of the Eq. (2.23) we obtain 

( (.Xp 5 Xp ) ) ( (Xp i Xp ) ) ■ 

Thus we find 

X~ \P P \ 2 = e 2 ^\a p \ 2 - 1. 

Now define 

[(XUX2)} := [(X1,X2)] (R) + [(XUX2)} {L) , 
which gives the contribution of the function \p to the total Klein-Gordon norm. Then we have 

K++ := [(Xp, Xp)} = ~e- np sinh 2 np (1 + X), 
K+- := [(Xp,X-p)} = ~e* p sinh 2 npa p f3 p , 

7T 

4 2 

K-+ := [(x-p,Xp)l = -e^sinh TTpa p P p , 

7T 



K— := l(x-p,X-p)} = -e*"smh 2 *p(l + X), 

7T 



and 



(2.25) 
(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 
(2.31) 
(2.32) 



(2.33) 



(2.34) 



To find the orthonormalized mode functions, we must find a matrix M„„i which diagonalizes the matrix K a 
(a, a 1 — +, — ).[] It can be chosen as 



K++ - A+ 



\ 



^\K + ^ + (K ++ - \ + Y v/|X + _| 2 + (K ++ A + )2 
-A-+ A++ - A_ 

V y/\K + -\* + (K ++ - \-)* v /|^+-l 2 + (^++-^) 2 / 



(2.35) 



Tnstead of allowing p to have the range — oo < p < oo, we restrict it in the range < p < oo but doubling the degrees of 
freedom with a = ±1. 



G 



where A± are the eigenvalues of K GIJ i, 

A± = \ ((K++ + X—) ± y/(K ++ - K^Y + 4|tf + _P ) . (2.36) 

Note that the matrix K aa -i is hermitian with positive eigenvalues, as it should be. Then the orthonormalized mode 
functions are given by 

V pa = 1 {M a+ Xp + M CT _ x~ P ) • (2.37) 



B. de Sitter super-curvature mode 

Next we consider the discrete mode that arises due to the spacetime structure of de Sitter space. As mentioned above 
we consider the wall perturbation mode seperately later. For convenience, we set p = iA and look for a possible value of 
A. When there is no bubble wall, it has been shown that there arises a mode with A = v' (y 1 = v/9/4 — to 2 / H 2 — 1/2) 
when v' > (to 2 < 2H 2 ) p2[ |. In the present case, because the wall is present and the mass changes from M to to, 
such a mode may or may not exist, depending on the function A4 2 even if m 2 < 2H 2 . 

In the region C, we consider the solution of the form, 

1 ^1-1/2 sinn r c ) 

u M m{rc,tc, Si) = -^r^ Xa(*c) . = Y lm (fl), (2.38) 

aE{tc) ^/icoshrc 

where we set the form of xa as 

= (a A P^(smtc-iO)+PAP-, A (smt c -iO), (t c -> -tt/2), 
XA lP i ; A (sint c - l O)+7 A pA( s i n t c „ i o), (t c -. tt/2), 

near the both boundaries of the region C. The regularity condition at tc = —it/2 is that the mode function should be 
less singular than (tc + 7T /2)~ 1 , hence xa ~^ 0. Using the asymptotic behavior of P£{z) at z — > — 1 p5[|, this requires 
the ratio of /3a to a A to be 

)1i = ^r(l+A + ^)r(A- A ')a A . (2.40) 

Similarly, the regularity condition at tc = vr/2 demands that 7a should vanish: 

7A = 0. (2.41) 



Then we need to solve the Eq. (2.18) with the above boundary conditions in both ends. Thus the problem is to find 
the eigenvalue A and the corresponding normalized eigenmode. This eigenvalue problem will be solved for a simple 
model in section 4. 

Now let us consider the normalization of the super-curvature mode. The Klein-Gordon inner products of the mode 
functions u\i m axe calculated in the region C as 



(uAlm {rc ,tc,ty, UM'm' {rc ,tc,ty) 



d UM'm' d UAlr, 



= -icosh r c J d£a E j dfl |u A ; TO — -^-um-t, 

X Sll'Smm' j (2.42) 



r(A + z + i)r(-A + z + i) 

where 



Na := I dC\xA\ 2 - (2.43) 



Thus, the normalized mode functions v\i m are given by 
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/r(A + ; + i)r(-A + ^ + i) 

VAlm = W — j^T" 1 WA/m . (2.44) 

In particular, by analytic continuation to the region R (inside the bubble), we obtain the de Sitter super-curvature 
mode functions in the open universe as 

VAim(t R , r R , SI) = -j= g- l sinh ^ yUm{r R , SI) , (2.45) 



where 



/r(A + 1 + i)r(-A + i + i) -P A -i/2 2 ( cosh m) v ... 

yAlm-=\ — -7T 1 , ■ , Y lm (Sl). (2.46) 

2 v / sinhr\R 



C. Wall fluctuation mode 



Recently, several authors have discussed the effect of the bubble wall fluctuations (T^JIJ,|T^]. The wall fluctuation 
mode appears when we consider the quantum flucuations of the tunneling field itself, which is the case of a single field 
model. As mentioned before, M 2 should be regarded as V"{(Tb) in this case. 

Since the wall fluctuation mode is one of the discrete modes, the formula obtained in the previous subsection is 
applicable. However, besides the fact that the origin o f the discrete spectrum is different from the de Sitter super- 
curvature mode, the corresponding solution of the Eq. ( 2.18| ), XWt can be written down formally in this special case. 
Hence we consider this mode separtely in this subsection. 

It is well known that the time derivative of the bounce solution satisfies the equation o f fluc tuations, related to the 
zero mode problem ]2q]. In fact, we can derive the following equation directly from Eq. (2.10), 



lP_d_ 

a% dtc 



dtr 



- V"(a B ) 



d<r B 
dt c 



0. 



(2.47) 



This is just the equation of fluctuations, i.e., eq.( 2.11 ) with the eigenvalue of — L 2 given by 1 + p 2 = —3, or A 2 = 



—p = 4. Then if they have finite Klein-Gordon norms, they will contribute to the quantum fluctuations inside the 
bubble. 

Therefore let us consider the mode functions 



da i 



Uw,lr, 



-l-l/2 
3/2 



Calculating the Klein-Gordon norms, we obtain 



(u\V,lm , UW,l'm< 



(i sinh rc) 



d(H 1 tc) y/i cosh rc 



rr(si). 



2AA 



r(/ + 3)r(i-i) 



Sll'Sri 



(2.48) 



(2.49) 



where 



w 



H 



+ 7T/2 



dtc a-E 



-x/2 



da i 



dtc 



+ 00 



da B 



(2.50) 



Note that Nw has the dimension of (mass) 2 . As one can readily sec, the Klein-Gordon norms vanish for / = and 
1. One may think that they give a divergent contribution to the fluctuations. However, they simply represent the 
temporal and spatial translations of the origin of the bounce solution, i.e., they are the zero modes. Hence we do not 
have to take these modes into account p6| . On the other hand, the modes with / > 2 do have finite Klein-Gordon 
norms and they contribute to the fluctuations pj| . We note that in the exact thin- wall limit, the amplitude of these 
modes are non-zero only on the bubble wall and they do not contribute to the fluctuations in the open universe (region 
R). It is crucially important that a B is not exactly zero off the wall, however small it may be. Then the normalized 
mode functions are given by 



VW,lr, 



'r(; + 3)r(; - 1) 



2M 



■ Uw,lm 



W 



(2.51) 
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Analytically continuing the mode functions into the region R, we have in the open universe, 



Vw,lm(r R , t R , Q) = —^=-^^—^y A=2Mn {r R , fl), (2.52) 



where J^AZm has been defined in Eq. ( 2.46| ). An important property of this mode is that the second rank tensor 



constructed from the spatial harmonic function happens to be transverse-traceless flll|| . Furthermore, the gauge- 
invariant density perturbation A, which represents the density perturbation on the comoving hypersurface, vanishes 
identically. Thus it may be viewed as a kind of gravitational wave mode [ Q . However, in order to treat all the modes 
on an equal footing, we do not take that view here. 

To summarize this section, we have discussed the three distinct kinds of mode functions which describe quantum 
fluctuations in open inflation models. The first one is the usual modes having continuous spectrum. The other two 
are modes with discrete spectra; the de Sitter super-curvature mode and the wall fluctuation mode. These quantum 
fluctuations give rise to cosmological metric perturbations which will be reflected in observational quantities. We will 
discuss this aspect of the quantum fluctuations in the next section. 



III. INITIAL FLUCTUATIONS IN OPEN UNIVERSE AND LARGE ANGLE CMB ANISOTROPIES 

Once we have the mode functions which describe the quantum state of the inflaton field inside the bubble, we can 
investigate the evolution of resulting cosmological perturbations in the open inflationary stage and the subsequent 
stage of the open Friedmann universe, which can be compared with observational data. In this section, we first 
consider the initial condition of the cosmological perturbations generated during the open inflationary stage, based on 
the results obtained in the previous section. Since the inflaton is almost massless during the open inflationary stage 
inside the bubble, we set m 2 = 0. Then we relate it to the spectrum of large angular scale CMB anisotropics. Below, 
we work in the region R which describes the open universe inside the bubble. We write the line element as 

ds 2 = -dt 2 + a 2 (t)^ i jdx q, dx j 

= a 2 (r]) -dw 2 +'f i jdx i dx j , (3.1) 

where a is the scale factor, t is the cosmological proper timej^] r\ is the conformal time, and 

j ij dx i dx j = dr 2 + sinh 2 rdfl 2 . (3.2) 
At the open inflationary stage, the scale factor is assumed to be given by 

s'mh Ht 1 

° = H = Fsinh(-77) ' (3 - 3) 

where H is sufficiently slowly varying in time. 

Just after the bubble is formed, or in the spacetime region close to the boundary light cone of the region R, the 
universe is dominated by the curvature term, or the matter energy density can be neglected. Hence the metric 
perturbations induced by the quantum fluctuations of the scalar field are negligible. On the other hand, there exists a 
time slicing on which the effect of metric perturbations is minimal, called the flat hypersurfaces. This implies that the 
fluctuations described by the mode functions we have obtained in the previous section may be regarded as those on 
the flat hypersurfaces. It is easy to show that the curvature perturbation on the comoving hypersurface 1Z C is related 
to the scalar field perturbation tp on the flat hypersurface as 1Z C = —d/(a^>B) t f p7| . Assuming the perturbation is 
adiabatic, 1Z C completely determines the metric and matter perturbations of scalar type (i.e., those which can be 
expanded in terms of scalar harmonics on the 3-space) . An important property of 1Z C is that it remains constant on 
superhorizon scales as long as the perturbation is adiabatic, irrespective of the background spatial curvature of the 
universe (see Appendix A). Hence what we need to evaluate is the amplitude of 1Z C for each mode when the mode 
passes the horizon scale. Assuming 4>b and H are sufficiently slowly varying at the inflationary stage, we may then 
evaluate 1Z C at the limit i] — > (a — > oo). Thus we have 

K c = —^-<p; r,^0. (3.4) 

<PB 



"•"Up to now we have used the non-dimensional time normalized by the Hubble parameter. But from now on, we recover the 
dimension of time in t. 
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The spectrum of the curvature perturbation on the comoving hypersurfaces is given as follows. Expanding the 
curvature perturbation by modes 



poo 

Tic = / dp U P aY p i m (r, Q) 

alm J ° 



(3.5) 



/■(. 



I 1 2 i 1 2 I 1 2 

the power spectrum of the curvature perturbation is given by \R- V \ '■= ^2 a \7Z- p(J \ for continuous modes, by 7?-a for 

I 1 2 

de Sitter super-curvature mode, and by \R-w\ for wall fluctuation mode. 



First consider the continuous modes. From Eqs. (2.37) and (3.4), we find after some algebra, 

coth up 



where 



Y 



H 2S 2 



b B J 2p(l+p 2 ) 



r(2 - ip) e-*P/3 p 



(i-n 



r(2 + ip) 2cosh7rpaj; 



c.c. 



(3.6) 



(3.7) 



It may be worthwhile to compare this result with the spectra of the Bunch-Davies vacuum state §Q and the 
conformal vacuum state fl27|,E8| . In these case we have 





(R 2 


I 2 - 
\'^p\bd ~ 


\^b~ 








(H 2 


1 P \ c 


\Jb~ 







coth 7rp 
L! J 2p{\+p 2 

2 1 



B ) 2p(l+p 2 



(3.8) 
(3.9) 



respectively. The result of our model (3.6) differs from that of the Bunch-Davies one in the last factor (1 — Y), which 
represents the effect of the mass difference outside the wall. However, as Y — > for p — > oo, the spectra in all the 

three cases behave as \7ip\ oc 1/p 3 for p> 1, which shows the perturbations have the Harrison-Zel'dovich spectrum 
on scales smaller than the curvature scale, irrespective of the choice of vacuum. The difference appears on and above 
the curvature scale, and it may become important when considering large angular scale CM B anisotropics. 
As for the de Sitter super-curvature mode, the spectral amplitude is given from Eq. (2.45) as 



1 



7V A r(2 + A) 2 ' 

We note that, in the Bunch-Davies vacuum limit, we have A — > 1 and N\ 



1/2, which gives 



BD 



(3.10) 



(3.11) 



The amplitude of the wall fluctuation mode is given from Eq. ( 2.52 ). We identify a B with <ps and note that 



daB/dtji — H &b- Then 



\1Z w\ 



H 2 



(3.12) 



One sees it is quite different from those of the other modes; it contains no information of how the field evolves at 
the open inflationary stage but is completely determined by the part of the potential which governs the false vacuum 
decay. 

Now we consider the CMB anisotropics predicted in our open inflation models using the initial conditions obtained 
in the above. The power spectra of CMB anisotropics in open models have been well investigated in the case when the 
scalar field is in the conformal vacuum [^-Q or in the Bunch-Davies vacuum [jsljisll. Since we expect the difference 
between the present models and the previously investigated models arises only on large angular scales, we focus on 
low mutipoles (I < 20) of the spectrum. Concerning the wall fluctuation mode, as we have mentioned, it has a very 
interesting property that the perturbation in the spatial curvature induced by it is transverse-traceless |ll| . Hence it 
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may be regarded as a tensor mode perturbation. Garriga have investigated the evolution of this mode in the open 
inflationary universe from this point of view H]. However, in this paper we treat it as a usual scalar perturbation 
except for the fact that p 2 = —4. For completeness, a proof that both approaches are equivalent is given in Appendix 
B. 

From the perturbed Einstein equations, the evolution equation of 1Z C becomes 



K" + 2—n' - KTZ C = S, 



(3.13) 



where the prime denotes the conformal time derivative, K = — 1 for an open universe and the source term S can 
be neglected if the universe is matter-dominated (see Appendix A). Since 1Z C — const, at the very early stage of the 
hot Friedmann universe, the relevant solution of it is the so-called growing mode solution. Here we consider 1Z C on 
sufficiently large scales that it re-enters the horizon after the universe becomes matter-dominated. In this case, the 
line element is described by Eq. (3.1) with the scale factor a(jj) oc cosh 77 — 1 (77 > 0). Then we obtain 



K r = He 



3(77 sinh ?y — 2 cosh n + 2) 
(cosh 77 — l) 2 



(3.14) 



where 1Z c j is t he ini tial a mplitu de determined at the inflationary stage of the universe, the spectrum of which is given 
by Eqs. @, jOc|) and 

Given the solution (3.14), we can compute the temperature fluctuations in the open universe in the gauge-invariant 
formalism plj . For this purpose, we consider the gauge- invariant variables \& and <£> which describe the gravitational 
potential perturbation and the curvature perturbation, respectively, on Newtonian slicing |H| . We have 'J + $ = 
in the present case. Then <I> is related to 1Z C as (see Appendix A) 



-K' 



(3.15) 



which gives 



$(ri)=1l Cti F(ri), 



(3.16) 



where 



3 (sinh 2 77 — 377 sinh 77 + 4 cosh 77 — 4) 
(cosh?7 — l) 3 



(3.17) 



Note that we have $ = t1Z c in the limit 77 — > 0. The potential perturbation gives rise to CMB anisotropies on large 



angular scales. Writing the temperature autocorrelation function in the form, 

C(0) = ^-5^(2Z + l)QP,(cosd), 



(3.18) 



the multipole moment C; is given by a sum of contributions from the continuous modes and the two super-curvature 
modes, 



Ci = cf ] + C\ a > + C\ 



(S) 



-,(W) 



(3.19) 



Assuming that the conformal time of the last scattering surface is sufficiently small, wlss ^ 1) but the universe is 
well matter-dominated by that time, we obtain 1 3 1 1 



C, (s) = \TZ A \ 



dp\TZ p 



F(r] LS s)fpi(va ~ Vlss) + 2 



dF(rf 



a 



(W) 



\TZw\ 



■^F(r] LS s) f\i {vo ~ Vlss) + 2 



:F(r) LS s)f2i(vo ~ Vlss) 



VLSS 

no at? 1^1 



drj 



— fpiivo ~ v')dv' 



dFW - ; 

' fAiiVo - V )dv 



Vlss 

710 dF ^) f i / U ; 

Vlss 



drj' 



(3.20) 
(3.21) 
(3.22) 



where 770 is the conformal time at present, f p i is the function given in Eq. J2.15D and 
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j a;= / r(A + j + l)r(-A + 1 + 1) P A U% 2 (<™kr R ) 

V 2 \J sinh rji 

Here a word of caution is appropriate. The redshift of the matter-radiation equal time is z eq w 4.2 x 10 4 i7o^ 2 , where 
ft is the Hubble constant normalized by 100 km/sec/Mpc and h < I from observations. In addition, coshr/Lss = 
1 + (Oq 1 — 1)/(1 + 2lss) where zlss is the redshift of the last scattering surface (~ 1000). Therefore the above 
formula for Ci becomes inaccurate for Qq -C 0.1. We should also note that we have put t]lss = for simplicity in 
actual calculations. Hence our results presented in the next section may have errors of ~ 10% for f2o ~ 0.1. 

Before closing this section, we mention that the formalism developed in the previous two sections can be extended 
to more general cases. Namely, as long as the regions R and L are described by de Sitter space, it is not necessary to 
have the same Hubble parameter for these two regions. Let Hr and Hl denote the Hubble parameters in the region 
R and L, respectively. Then it is easily recognized that all H 's appearing in the formulas in the present section should 
be identified with Hr. Further, in this case, as the scale factor is different from that of de Sitter space in the regio n 



C, but the formulas written in terms of the conformal coor dinate £ do not contain H explicitly except for Eq. (2.18) 
what we need to do is to replace the potential in Eq. (ETq) by 



U(0^-^f~l + a%M 2 , (3.24) 
a E d£ 2 



where £ is now defined in such a way that the metric in the region C takes the form (2.17) with Oe(£) being a general 
function of £. 

IV. SPECTRUM OF THE CMB ANISOTROPIES PREDICTED BY SOME MODELS 

In this section, we apply our formulas obtained in sections 2 and 3 to simple models of the open inflation scenario. 



A. Single field model 



We first consider a model with a single scalar field with the potential as shown in Fig. 1(a). The tunneling 
field becomes the inflaton field after the tunneling. Thus the mass square of the fluctuation field ip is given by 
A4 2 (t) — V"{<jb)- Rs typical shape is rather complicated as shown in Fig. 3. 

Because of the complexity of M 2 (t) 1 the analysis of this model is difficult in general. However, the mass of the 
inflaton in the false vacuum M must be large compared with H in order that the tunneling is not dominated by 
the Hawking-Moss instanton (33) as explained in Ref. 0. Furthermore, as there are many subtleties^] in models with 
M 2 ~ H 2 , we leave it as a future problem. Here we assume M 2 ^> H 2 . 

Then the potential barrier between the two asymptotic regions at £ — ► ±00 is so high that the coherence between 
regions R and L are exponentially suppressed as long as we consider the modes with p 2 <C M 2 /H 2 , which is the case 
of our present interest. So one of the two modes corresponding to the same p can be set to vanish at £ — > — 00 and 



the other at £ — > +00. Then automatically they are orthogonal to each other. Since the field equation (2.18) is real 
in the region C, we can choose the mode which vanish in the region L to be real in the region C. Then analytically 
continucing it to the region R, the normalized mode which does not vanish in the region R can be written as 

^ - ( ssfe (g±i)%.-*-*(« + ») (gii)-*) . M , 

( 0, (in L) 

where we set m — and 5 P is a real constant which depends on the detail of the potential. Using this mode function, 
the amplitude of 1Z P is evaluated as 

It? i 2 _ ( H2 \ 2 cosh 7rp + cos 2<5 P 

1 H ~\&b) 2p(l+p 2 )smh7rp { ' 



§ For example, the formalism developed in |17J which we used as a given result in the present paper is not applicable to the 
case in which the wall is spreading so broadly that we cannot approximate o"s(r) by of for r < 0. When M 2 ~ H 2 , the wall 
necessarily becomes broad. 
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Comparing this result with the spectra for the Bunch-Davies vacuum (3.8) and the conformal vacuum (3.9), we find 

■ 1 2 I 1 2 

that the difference between the present spectrum and \R- P \ BD is no more greater than the difference between 



and 1 72.p | ^ . As was clarified in the effect of the difference between the Bunch-Davies vacuum and the conformal 

( c) 

vacuum on C) is always negligiblly small independent of £Iq. Thus as far as the continuous modes are concerned, 
we do not have to perform further calculations. 

Next we consider the discrete mode. As mentioned above, the wall fluctuation mode is always present when we 
consider the flucuations of the tunneling field itself. Then the question is whether there exists another discrete mode 
like the de Sitter super-curvature mode. The answer is no as long as M 2 > Hi. The reason why is explained in the 
Appendix C. In the general case, we do not have an answer to this question,[J but we do not discuss the de Sitter 
super-curvature mode any further. We focus on the wall fluctuation mode below. 

The wall thickness is roughly evaluated by the inverse of the curvature scale of the potential as M . Since we are 
now assuming that M 2 ^> H 2 , it will not be too crude to adopt the thin- wall approximation. Then one can evaluate 
Afw as 



w 



H J dTa E (r)(^) 2 ^RS i; 

S i:=H Jdr(^) 2 , (4.3) 

where S± is the tension of the wall. If one uses the reduced Euclidean action written in terms of R, S\, and the potential 
energy density inside the bubble Vr (true vacuum) and outside the bubble Vl (false vacuum), R is determined in 
terms of the other parameters by minimizing the action as |37| , |l5|| 

R= - 3Sl (4.4) 

y/(AV + QttGS 2 ) 2 + 24ttGV r S 2 

where AV — Vl — Vr. Then the final result for the amplitude of the wall fluctuation mode is 



\U W \ 2 = (AV^ + 6nGS 2 ) 2 + 2A7:GVrS 2 



H 



8ttGS 2 



(H 2 -H 2 + {AttGS^Y + {^GS^r) 2 , (4.5) 



where Hl and Hr are the Hubble parameters in the regions L and R, respectively. Note that although we assumed 
(Hl — Hr)/ Hr <C 1, the difference Hl — Hr may not be negligible in the determination of the amplitude of the wall 
fluctuation mode. 

Now we consider the CMB power spectrum due to the wall fluctuation mode. Fig. 4 shows the power spectrum of 

CMB anisotropies 1(1 + \)C[ W) for various values of Slo ■ The curves are normalized by H 2 / RS\. As discussed before, 
the amplitude of the CMB anisotropies due to the wall fluctuation mode are determined by the nature of the bubble 
wall. This gives an independent constraint on open inflation models. In particular, the spectrum of CMB anisotropies 
rises sharply towards lower multipoles for open models with Oo ^ 0.2. Hence if this mode dominates over the other 
mode contributions, such single field models will be severely excluded by the COBE data pq ]. This issue has been 
recently investigated by Garriga |l4| and Garcia-Bellido 

To carry on a further analysis, we have to specify the model in more details. Let us parametrize the potential V(a) 
by M 2 , V L , V b and AV as shown in Fig. 1(a). We assume V(a) < Af* j5 M < M pt where M pl is the Planck mass, 
AV/Vl -C 1 and the potential has the unique curvature scale M 2 in the region C. Then, as the wall thickness is 
essentially given by M _1 , we have 

Further we assume V(a) ~ XM* t (a / M p i) 2n (n = 1, 2, • • •) on the right of the barrier, where A <C 1 as in the chaotic 
inflation. 



**We note an issue related to this problem. One can show that the non-existence of another discrete mode in the range 
-4 < p 2 < is a sufficient condition for the uniqueness of the negative eigenvalue mode in the one-loop order calculation of 
the tunneling rate in the path integral approach |34|-|3q] . But the uniqueness has not been proved for the tunneling on the de 
Sitter background as far as we know. 



13 



There are three requirements to be satisfied. 1) The tunneling must be dominated not by the Hawking-Moss 
instanton but by the Coleman-De Luccia instanton. This requires that M 2 > 4H 2 . 2) The tunneling rate T must be 
suppressed enough in order to avoid the fluctuations caused by the bubble collisions. This implies that [BIS 



- ln ( r ^ 4 )-^» L M 

3) Finally, the wall fluctuation mode must not dominate 1Z C . Noting that the contribution of the continuous modes 
to its power in the logarithmic interval of p at p ^S> 1 is (1Z 2 ) := liirip^oo \lZ p \ 2 p 3 / (2ir 2 ), this requires 



p- 

2 x 2 



l^l 2 = (&) 



(4.8) 



We assume the requirement 1) is satisfied and investigate the conditions on the potential parameters derived from 
the requirements 2) and 3). 

From the slow roll equation of motion, we have \H^/&b\ = 8itHrVr/\V' (cb)| = $>n(H r/M p i, where £ > 1. 
Therefore 



\Kw\ 2 7T \( Ay \ 2 / H R 



(R 2 ) 8C 2 V V GttGS 2 ) \2nGSi 
Since the left hand side of this must be smaller than unity, we should have 



+ UzT^r • (4-9) 



Ar <C 2 and ^V<C 2 . (4.10) 



QttGS 2 ~ s 2ttGSi 

Using the fact Si ~ Vh/M and AV -C Vl, we find the first inequality is automatically satisfied if the second one is, 
which is re-expressed as 

V b > HrMM 2 [ ^ 4 (H R \M 

V L ~ 2irC 2 V L 3C 2 \H L ) H L - ( ' ' 

If this condition is satisfied, we also find the requirement 2) is fulfilled. Thus if the potential barrier is high enough, 
approximately of 0(M/Hl) x Vl, the wall fluctuation mode will become harmless. Apart from the intrinsic unnat- 
uralness of a single field model, it is not difficult to construct models which satisfy the above constraint. Hence we 
conclude that a single field model remains still viable. 

B. Two-field model 

Next we consider a two-field model in which the tunneling field a and the inflaton in the nucleated bubble (f> are 
different. The supernatural inflation model proposed by Linde and Mezhlumian Q is included in this category. We 
consider the following situation. Before tunneling, <p is at the minimum of the potential cf>L with mass M. During 
the tunneling, the potential of <j> changes to an almost flat but slightly declined one, i.e., Vi ^ and V£ ~ 0. Thus 
4>l is no longer the minimum of the potential in the nucleated bubble. Therefore <j> begins the slow rolling to the new 
minimum of the potential. 

In the following calculation, we assume that the gradient of the potential is so small that we can neglect its effect 
in estimating the amplitude of the fluctuations tp, though the condition under which this neglection is justified is not 
clear. Then our present formalism can be applied. For simplicity, we also assume that the thin wall approximation 
is valid. We denote the radius of bubble wall by R = ae(ro) = H^ 1 cos To. Then the bubble wall trajectory in the 
Lorentzian region is described by the hypersurface tc = T o = const, in the region C. 

Under these assumptions, we can evaluate the fluctuation spectrum inside the bubble by applying the formulas 
obtained in section 3 with M 2 given by 

M 2 -{ 0; (to < *c < tt/2), 



The form of the potential U(£) in Eq. (2.18) in this case is illustrated in Fig. 5. Thus the parameters of the present 
open inflationary model are the mass of the inflaton field in the false vacuum M, the wall radius R and the Hubble 
parameter H. 
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First we consider the cont inuous spectrum. As the mass is constant both outside and inside the wall, we find that 
the solution of the Eq. (|2.1§| ) becomes 



4 R) (t c ) 
4 L) (t L ) 



^(sinic 



iO); {t <t c < tt/2) 
iO); (-tt/2 < t c < T Q ). 



(4.13) 



The coefficients a p and j3 p are determined by the junction condition at the wall in the region C: 



dxi R) 



( X% 



Xp 



dx ( - L) P d Xp L) 



V di J 



(4.14) 



Then a p and /3 P given by Eqs. (2.28) can be now easily evaluated on the wall. 

I 1 2 

We have numerically evaluated the spectrum of the continuous modes. Figs. 6 show the power spectra of \R- P \ 
for various values of M and R, normalized by the Bunch-Davies vacuum spectrum. Fig. 6(a) shows the case when 
the mass outside the bubble is M/H = 2. The lines are the power spectra for the wall radii HR — 0.1, 0.5, 0.7, and 
0.9. One sees the spectra for all the wall radii are almost the same and they coincide with that of the Bunch-Davies 
vacuum except for the small range of very small p. On the other hand, Fig. 6(b) is the case M/H = 10 with the 
lines showing the power spectra for HR = 0.1, 0.5, 0.7, and 0.9. For M/H 3> 1, we can apply the discussion given 
in section 4A, around Eq. ([T^). There we have seen that the spectrum does not differ much from the case of the 
Bunch-Davies vacuum. Thus the increase in the amplitude at p < 1 will saturate as M/H — > oo. In both Figs. 6(a) 
and (b), all the lines rapidly approaches unity at p 1, which corresponds to the Harrison-Zel'dovich spectrum, in 
accordance with our expectation. 

Next we turn to the discrete spectrum. Again the mode funct ion is readily solved by means of the associated 
Legendre funciton. Thus eigenvalue probem reduces to Eq. ( p. 40 ) with a a and /3a given in terms of the junction 
condition at tc = tq as 



((V>^ A )) 



OA = 



0A 



*c=n> 



(4.15) 



It is instructive to show an analytically solvable example of this eigenvalue problem. Let us consider the case 
M 2 = 2H 2 (/it' = 0) and m 2 = (y' = 1). In this case, the associated Legendre functions are expressed in terms of 
elementary functions. 



Pt\z) 



1 



1 



T(1=f A) \z - 1 



z =f A 



±A/2 



±A/2 



T(2t A) \z- 1 

Also, Eq. ( 2.40 ) reduces simply to /3a = 0, hence from Eq. ( [4.15 ), 



((P-\P A ))=Q. 



This is easily solved to give 



A : 



-Xq + y/2 - X% 



(4.16) 



(4.17) 



(4.18) 



where xq = sin tq = \/T— (HR) 2 . This result shows that A becomes larger as the wall radius increases. 

For general M 2 , we need numerical evaluation. Since < A < 1 in the exact de Sitter case with the upper limit 
attained when m 2 — > 0, we expect any introduction of finite mass outside the wall will reduce A to a value less than 
unity. Hence we look for an eigenvalue in the range < A < 1. Fig. 7 shows A as a function of HR for several 
values of M 2 when m 2 = 0. One sees A becomes smaller as M 2 becomes larger or as the wall radius becomes smaller. 
But the de Sitter super-curvature mode exists for any wall radius as long as M 2 < 2H 2 . On the other hand, for 
M 2 > 2H 2 , the de Sitter super-curvature mode ceases to exist when the wall radius becomes smaller than a critical 
value. Fig. 8 shows the critical line on the (M/H, i/i?)-plane on which the super-curvature mode disappears. The de 
Sitter super-curvature mode exists below the line. 
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This property of the super-curvature mode can be understood in analogy with the quantum mechanics described 
before. As the wall radius decreases, the position of the bubble wall £o moves right in Fig. 5 of the potential f (£). Also, 
as the mass M increases, the potential barrier outside the bubble wall £ < £o becomes higher. These make it difficult 
to form a bound state. We note that U(£) has the deepest valley when the mass M is zero. This situation corresponds 
to the case when we assume the Bunch-Davies vacuum state inside the bubble and when the super-curvature mode 
contributes extremely [^). 

Now we show the results of numerical calculations of the CMB anisotropy power spectrum 1(1 + l)Ci in Fig. 9 for 
the f^o = 0.1 universe with M 2 /H 2 = 2 and the wall radii HR = 0.3, 0.5, 0.7, and 0.9. The curves are normalized by 
(3Jf 2 /5(/)B) 2 . Here both continuous and discrete modes are included. To compare the results with the previous ones 
p7|-p0|,p|Jl^1 , we plotted two dashed lines in the figure. The top dashed curve is the result when the scalar field is in 
the Bunch-Davies vacuum |||l3), and the bottom dashed curve is the one in the conformal vacuum 27- 30|. All of the 
curves lie between these two curves. We see the amplitude becomes large as the wall radius increases, approaching 
that in the case of the Bunch-Davies vacuum. As described before, the contribution of the super-curvature mode is 
most for the Bunch-Davies vacuum and is least for the conformal vacuum. When the wall radius becomes large and 
the mass becomes small, the contribution from the super-curvature mode becomes large, which explains the behavior 
of the CMB power spectra in Fig .9. We have also calculated the CMB spectra for the f2o = 0.3 universe. The results 

(c) 

have turned out to be almost independent of the model parameters. In fact, if only Cy is taken into account, the 

difference between the Bunch-Davies vacuum and the conformal vacuum is only a few percent even for Oq = 0.1 [0. 

(S) . . . . 

Thus the differences are dominantly due to the de Sitter super-curvature mode (7, , but its contribution rapidly 
becomes negligible as one increases SIq- 



V. CONCLUSION 



We have investigated in detail the quantum fluctuations and the resultant CMB anisotropies in two simple models 
of one-bubble open inflationary scenario. Before going into the analysis of the specific models, we have extended the 
previously developed formalism to investigate the quantum state inside the nucleated bubble, and derived formulas for 
the mode functions which take account of the effect of the tunneling described by the Coleman-De Luccia instanton. In 
the main stream, we have assumed the potential energy difference between the false and true vacua is small enough so 
that spacetimc both inside and outside the bubble can be described by de Sitter space with a single Hubble parameter 
H, but the results can be easily extended to more general cases. 

A complete description of a quantum state requires a complete set of mode functions which are normalizable on 
a Cauchy surface of the whole spacetime. This brings about new sets of fluctuation modes which had not been 
considered in most of the previous analyses of an open inflationary universe. In addition to the usual modes with 
continuous spectrum, there are modes with discrete spectrum corresponding to the fluctuations on super-curvature 
scales. For both the continuous and discrete modes, we have obtained general formulas of the spectrum of the CMB 
anisotropies on large angular scales which result from the quantum fluctuations. 

After these preparations, we have considered two simple classes of one-bubble open inflationary models. One is a 
single field model with the potential as illustrated in Fig. 1(a), assuming that the mass square in the false vacuum, 
M 2 , is much larger than H 2 . The other is a two-field model in which the false vacuum decay is mediated by a scalar 
field different from the inflaton inside the bubble and the inflaton is massive in the false vacuum through the coupling 
with the tunneling field. 

In the case of a single field model with M 2 ^> H 2 , there is one discrete mode which represents fluctuations of the 
bubble wall. We have shown that there is no other discrete mode in this model. Thus we have considered the CMB 
anisotropies due to the continous modes and the wall fluctuation mode. As far as the contribution of the continuous 
modes to the CMB anisotropies is concerned, we have found that it is approximately the same as in the case of the 
Bunch-Davies vacuum or the conformal vacuum. Hence we have focused on the wall fluctuation mode. The curvature 
perturbation due to the wall fluctuation mode is notable. It is always present in this model and its amplitude is 
totally determined by the part of the potential which governs the tunneling but has nothing to do with the scalar 
field dynamics inside the nucleated bubble. Since its amplitude is independent of the amplitude of the curvature 
perturbation due to the continuous modes, its contribution to the CMB anisotropies gives rise to an independent 
constraint on the model, just as the gravitational wave perturbation does to the usual inflation models. Fortunately, 
the constraint turns out to be relatively weak. Thus it is possible to construct a model in which the effect of the wall 
fluctuation mode can be neglected. 

Here we make the following remark. The wall fluctuation mode has a peculiar property that curvature perturbations 
induced by it are transverse-traceless, hence it can be regarded as a mode of tensor- type perturbations. This means 
that the perturbation of the scalar field might couple with the gravitational wave perturbation for this mode when 
the degrees of freedom of gravity are fully taken into accout. But we have not incoorpolated them in the present 
analysis. Thus there remains a possibility that the final answer changes qualitatively. We hope to come back to this 
issue in near future. 

As for a two-field model, we have adopted the thin-wall approximation and assumed the mass of the inflaton changes 
like the step function across the wall. Thus the model can be parametrized by the mass in the false vacuum M, the 
Hubble radius H and the wall radius R. In this model, the origin of the discrete mode is different from that in the 
single field model. Since it originates from the spacetime structure of de Sitter space, we have called it the de Sitter 
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super-curvature mode. The existence of the de Sitter super-curvature mode depends on the mass M of the inflaton at 
the false vacuum and on the radius R of the bubble wall. It appears when HR is large and M/H is small. For models 
with M/H 1, the super-curvature mode disappears. Different from the wall fluctuation mode which appeared in 
the single field model, the amplitude of spatial curvature perturbations induced by this mode is determined by the 
scalar field dynamics at the open inflationary stage inside the bubble, as in the case of the continuous modes. 

We have then investigated the spectrum of CMB anisotropies on large angular scales in this model. Though the 
spectra of the curvature perturbations due to the continuous modes have different shapes on the curvature scale 
for different model parameters, we have found the resulting CMB power spectra do not significantly depend on the 
parameters for Slo 0-1. On the other hand, the parameter-dependence appears clearly in the contribution of the 
de Sitter super-curvature mode to the CMB spectrum. The effect of this mode is appreciable when the wall radius is 
large (HR ~ 1) and the mass is small (M/H <C 1). The effect is to raise the amplitudes of low multipoles at I < 10 
for models with Hq < 0.1. If this mode contributes significantly, the predicted CMB power spectra will contradict 



with COBE observations |13|. However, for an open universe of f2o ^ 0.3 as well as for models with small HR and 
large M/H, the effect of the super-curvature mode on the CMB anisotropies is practically negligible. 

We comment on implications of our results to the large scale structure formation. As all our open models predict 
the Harrison-Zel'dovich spectrum on small scales, the difference can appear only from the normalization of the density 
perturbations. However, the difference will be negligible if we adopt the normalization scheme in terms of the likelihood 
analysis using the COBE result (3^] . This is because all the models predict practically the same CMB power spectrum 
at I > 10. 

In summary, we have presented a detailed analysis of the quantum fluctuations in open universe in simple models of 
the one-bubble inflationary scenario and the resulting CMB anisotropies on large angular scales. We have found there 
exist ranges of model parameters which are consistent with CMB observations. Thus the one-bubble inflationary 
scenario is a viable one for explaining the large scale structure of the universe and it certainly deserves further 
study. There are of course many issues left to be clarified in future. For example, consideration of the continuous 
gravitational wave modes is definitely necessary. Analyses of more sophisticated models, like several other two-field 
models proposed by Linde and Mezhlumian are of particular interest. Inclusion of the degrees of freedom of the 
gravitational perturbation in the formalism from the beginning is a difficult issue but should be done in order to gain 
a more firm picture of the false vacuum decay and the subsequent evolution of the quantum state. 
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APPENDIX A: DERIVATION OF THE EVOLUTION EQUATION FOR TZ C 

Here we derive the evolution equation for TZ C . We follow the notation of Ref. [B2[ for the perturbation variables. 
Consider the general scalar harmonics on 3-space with curvature constant K, whicn satisfy 

(3) 

(A+fc 2 )F = 0. (Al) 
The corresponding vector and tensor harmonics are defined as 

Y * = ~l Y \i> Y » = ^ Y \ij + l^ Y (A2) 

where the vertical bar denotes the covariant derivative with respect to the 3-metric 7y . In the notation of the present 
paper, we have set K = — 1 and k 2 = p 2 + 1. But since we do not have to fix K in the discussion below, we leave it 
arbitrary and use k to denote the eigenvalue. Also we suppress the eigenvalue index k for notational simplicity. We 
expand all the perturbation variables in terms of these harmonics. Specifically, the metric is expressed as 

ds 2 = a 2 [-(1 + 2AY)drf - 2BY i dx i dr 1 + ((1 + 2H L )^ ij + 2H T Y lJ ) dx l dx 3 ] , (A3) 

and the energy-momentum tensor is expressed as 

T° = -p(l + SY), T\ = -(p + p)vY\ 

T\j =p ((l + ix L Y)5) + 7r T Y l 3 ) . (A4) 
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The background equations are 

&'\ 2 8itG n , a' , . , 

- +K=— pa 2 , p' + 3-(p+p) = 0, A5 

a J 6 a 

where the prime denotes the derivative with respect to r/. Then the Einstein equations for scalar-type perturbations 
are written as 

SG° = SnGpSY , 

8G° i = 8^G(p + p)(v-B)Y i , 

5G l , = 8wG(pir L Y 7ij + p^ T Y^ ) , (A6) 

where the explicit expressions for SG^v can be found in Appendix D of Ref. p2[ . 

For our discussion, we need the (0, i)-component and the traceless part of (i, j)-component of the above equations. 
They give, respectively, 

k—A - kW + Kcr g =4irG(p + p)a 2 (v- B) , (A7) 



k 2 (A + 11)- -^(a 2 ka„)' = -8TrGpa 2 7r T , (A8) 
or 



where 1Z = Hl + hHr and a g = k 1 H' T — B. Expressing Eq. ( A7) on the Newtonian hypersurface (defined by <r g = 0), 
for which we have A = ^, 1Z = v — B = V, we find 

$' = -ty - 4irG(p + p)a 2 ^ . (A9) 
a k 

Also, from Eq. (|A8|), we obtain 

a 2 

# + $ = -8ttG-^pII , (A10) 



where IT = tt^ is the gauge-invariant anisotropic stress perturbation. 

Using the contracted Bianchi identities, which give the energy-momentum conservation law T^ u - u = 0, one obtains 
the equations for the matter variables. Here we only need the p — i component of it, i.e., the equation for V. It is 
given by Eq. (4.76)' in Chapter II of ]3^| , which is 

1 ( aV \ T a c c ? A + r 2fc 2 -3if w _ , A11 , 

- T =^ + S; S:=^— n, (All) 

a \ k J 1 + w o k 1 + w 

where c 2 = p'/p', w — p/p, T is the gauge-invariant entropy perturbation, and A is the density perturbation on the 
comoving hypersurface which is related to that on the Newtonian hypersurface A s as 

A = A s + 3(1 + w)-^- . (A12) 
a k 

Now from the gauge transformation property of 1Z, the curvature perturbation 1Z C on the comoving hypersurface 
(defined by v — B = 0) is expressed in terms of $ and V as 

ftc = *--^. (A13) 
a k 



Then from Eqs. (A9) and (All), we find 



H'=-K^--—S. (A14) 
k a 



Taking the derivative of this equation and using Eqs. flAiol) and ( |All| ), we finally obtain 

a', / „ / a \ „ \ „ / a ' \ K 



K'Z + 2—Tl' c -Kn c = -l2l — \ +K\ S- I —S\ +^8irGpa 2 n. (A15) 
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We find the right hand side of this equation vanishes if the universe is matter-dominated (p — c 2 — 0). Furthermore, 
it is known that A = 0((ka/a') 2 ) on superhorizon scales |32|. From this fact, one can easily deduce 1Z C remains 
constant on super- horizon scales if T = II = 0. We also note that from Eqs. (A13) and (A14), $ can be expressed as 



K a a 

In particular, if S — as in the case of a matter-dominated universe, we have 

1 a! . 

<t> = n c ---n' c . 

K a 

With K = — 1, this is the relation we have used in the text. 



(A16) 



(A17) 



APPENDIX B: EQUAVALENCE OF TENSOR AND SCALAR TREATMENTS FOR P 2 = 4 

Here we show that the wall fluctuation mode can be treated either as a tensor-type perturbation or a scalar-type 
perturbation and prove the equivalence of the CMB anisotropy formulae in both approaches. 
In general, setting ST/T — O, the perturbed collisionless Boltzmann equation is written as 



D 



1 



— 6 = kAYri 1 - (K' - -ka g )Y - kagY^fj 



(Bl) 



where D/dX denotes the Lagrange derivative along the light ray with A being the conformal affine parameter, and 7; 
is the directional cosine of the photon propagation vector. Assu ming the universe is matter-dominated, one has A = 
on the comoving hypersurfacc v — B = 0. Also, expressing Eq. (A7) on the comoving hypersurface, one obtains 



K 



(B2) 



where {ka g ) c is ka g evaluated on the comoving hypersurface. Using these facts, Eq. (Bl) on the comoving hypersurface 
becomes 



n k 2 — "\K 

^X C = -^^—( k ^)cY - (ka g ) c Y l3 Yl J ■ 
In particular, for k 2 = 3K which corresponds to p 2 = k 2 — 1 = —4 for K = —1, one has 



(B3) 



(B4) 



Noting that in this case is transverse-traceless 11 , we see this is equivalent to the tensor- type perturbation with 
the metric perturbation Hij = 61Z c Yij. 

Now we go to the Newtonian slicing, in which a g = 0, by the transformation rj — > r/ + T. Then the gauge function 
T is determined from 



which gives 



0= {ka g ) c -k 2 T, 
(ka g ) c 



T 



k 2 



K 

K 



(B5) 



(B6) 



Let Q s denotes the temperature anisotropy on the Newtonian hypersurface. The gauge transformation affects only 
the monopole and dipole parts of the anisotropy. Then one finds <d s is related to 9 C as 



-TY + TYirf 



(B7) 



Inserting Eq. (B6) into this equation, we obtain 



[©J 110 = [Qs - - iKY - ^KYul 1 ] 

l L J»?iss l a K K 1 J 



no 

VLSS 



(B8) 
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where [• • ss denotes the difference between the quantity evaluated at t] and tjlss- 
On the other hand, the equation for 9 S is obtained from Eq. (Bl) by setting a g = as 

— 6 S = $l, l7 l - &Y = -tt($Y) - 2&Y , 



(B9) 



where we have used the fact $ + ^ = 0, which follows from Eq. ( A10 ). Combining Eqs. (BS) and (B9), we find 



I m 

' Vlss 



a' I , 1 „ n 

= <[, _ --K' C Y - ^KY\*l l ] 



a K 



Vo 

Vlss 



(BIO) 



VLSS 



Noting that we have 



-K' 



5 a c 5 

in the limit rj —* (with K = —1), the above equation reduces to 



[e c 



vlss 3 



2 / $'ydA + (monopole + dipole) . 

J VLSS 



(Bll) 



(B12) 



for r/Lss — * 0. This is just the formula to evaluate the CMB anisotropies due to scalar-type perturbations |3l]]. Thus 
we have proved the equivalence of the tensor and scalar approaches to the CMB anisotropies when p 2 = —4, i.e., for 
the wall fluctuation mode. 



APPENDIX C: ABSENCE OF THE DE SITTER SUPER-CURVATURE MODE IN THE SINGLE FIELD 

MODEL 



Here we show there exists no discrete mode other than the wall fluctuation mode in the single field model, provided 
the wall thickness, i.e., the inverse mass scale of the potential is much smaller than the Hubble radi us H ~ 1 as well 
as than the wall radius R. We consider the following simplified model in which the potential in Eq. ( [2. IS ) takes the 
form, 



Ui 



cosh Co 



u(0 = 



cosh £ 
U u 



cosh 2 £ 



(C < Co - e) ; (I) , 

(Co - e < C < Co) ; (II) , 

(Co < C < Co + e) ; (HI) , 

(£o + e < C < Co + e) ; (IV) , 



(CI) 



where Ui := (M 2 /H 2 - 2)/ cosh 2 Co (see Fig. 10). We assume U\, Uu, Um 3> 1, in accordance with the assumption 
M 2 /H 2 3> 1. We note that Co must be positive, since d is negative inside the wall. We consider the solution \v with 
p = iA, where < A < 2. Then we may approximate the potential in (I) by the constant U\. Further, for simplicity, 
we set U\ — Uu\. Then the solution takes the form, 



XA 



where 



Cie K («-«") (I), 

dl (e lfc («-«°) + &„e- 4fc «-«°)) (II), 

Cm (e K («-«°) + bme- K[ ^ o) ) (HI), 

Ci V e- A( ^ o) (tanhC + A) (IV), 



K = y/Ul + A 2 , k = y/Uu - A 2 , 



(C2) 



(C3) 



and C's and 6's are constants determined by the continuity conditions of xa an d dx\/dt;. In (I), we take only 
the growing solution (oc e K ^) because the contribution from the decaying solution (oc e~ K ^) will be exponentially 
suppressed. 

Solving the continuity condition of d(lnxA)/^C at C = Co — £ and C = Co, we find d(lnxA)/dC at C = Co + £ is 
calculated from the left of the point to be 



20 



1l> 



L ■- 



where 



sin(2# — ke)e K£ — sin ke e K£ 
sin(2# — ke)e Ke + sin ke e~ Ke ' 



tan 9 = n/k . 



On the other hand, the same quantity is calculated from the right of £o + e to be 

1 



-A 



coslU(£ + e)(tanh(£ + e) + A) 



Thus the continuity condition at £ = £o + £ reduces to 



(C4) 



(C5) 



(C6) 



(C7) 



This equation must be satisfied for A = 2 because we know the wall fluctuation mode always exists. This condition 
gives one constraint on the model parameters U\, U\\ and e. Therefore all of them cannot be chosen arbitrarily. This 
reflects the fact that M 2 = V"(ctb)- Suppose we fix U\ and Uu, or equivalently, k := k|a=2 an d k := k\^ =2 . Then s 
is determined from Eq. ( K7u with A = 2. In order to find the order of magnitude of s, we note that — 2 < tpR < — 3/2 
when A = 2, i.e, ipR = 0(1). Since k» 1, this implies the numerator in ip^ must be small, 



sin(2f? - ke) e Ke - sin ke e K£ < 1, 



(C8) 



where = 
A = 2 as 



= 2- Thus we find that e = 0(1/k). Using the above inequality, the denominator in tpi, is evaluated for 



2 sin ke e 



(C9) 



and is of order unity. 

Now we show that no other discrete mode exists. Suppose there existed another discrete mode. The next largest 
eigenvalue solution must have one node. Noting the constraints £o > and A > 0, we can easily see that there would 
be no node in (IV). Thus the solution must have a node in (II) or (III). As we decrease A from 2, because xa must 
vanish at £ = £o + e just before a node would first appear in (III), ipL would diverge for some value of A. However, the 

changes of k and n due to the variation of A would be at most of 0(l/k) and 0(1/ k), respectively, and accordingly 
the change of 9 would be also small. Hence there is no chance that the denominator in ipL would vanish, which is a 
contradiction. Thus we conclude that there is no discrete mode other than the wall fluctuation mode in this limiting 
case. 
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Figure Captions 



Figure 1. A schematic picture of the tunneling potential for the one-bubble open inflationary scenario, (a) for a single 
field model and (b) for a two-field model. 

Figure 2. Penrose diagram of the de Sitter space. The coordinates which cover the regions R, L and C are shown, 
respectively. 

Figure 3. A schematic picture of Ad 2 as a function of r (= tc) in the case of a single field model. 

Figure 4. The power spectrum of CMB anisotropies due to the wall fluctuation mode for various values of Hq. The 
lines are normalized by H 2 /(S\R). 



Figure 5. A schematic picture of f (£) in eq.(2.19) for a two-field model with thin-wall approximation. 

Figure 6. The power spectra of curvature perturbations due to the continuous modes. The real line, the long dashed 
line, the short dashed line and the dotted line are, respectively, for HR = 0.1, 0.5, 0.7 and 0.9, for (a) M/H = 2 and 
(b) M/H = 10. 

Figure 7. The eigenvalue A of the de Sitter super-curvature mode as a function of HR where R is the wall radius. 
The three lines, from top to bottom, show the cases for M 2 /H 2 = 5/4, 2 and 9/4, where M is the mass outside the 
bubble. 

Figure 8. The critical line on which the de Sitter super-curvature mode disappears on the (M/H, HR)--pl&ne. The 
super-curvature mode exists in the region below the line. 

Figure 9. The CMB anisotropy power spectra in the open universe with Slo = 0.1 predicted in two-field models. 
The lines are normalized by (3iJ 2 /5</>s) 2 . The top and the bottom curves are the results when the scalar field is in 
the Bunch-Davies vacuum and the conformal vacuum, respectively. For the other lines, the model parameters are 
M 2 /H 2 = 2 and HR = 0.9, 0.7, 0.5, and 0.3, from top to bottom. 

Figure 10. A simplified potential U(£) for a single field model. 
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